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Abstract 

High transverse momentum tt° -mesons have been measured with the HI detector at HERA 
in deep-inelastic ep scattering events at low Bjorken-x, down to x 4-10~ 5 . The mea- 
surement is performed in a region of small angles with respect to the proton remnant in the 
laboratory frame of reference, namely the forward region, and corresponds to central rapid- 
ity in the centre of mass system of the virtual photon and proton. This region is expected to 
be particularly sensitive to QCD effects in hadronic final states. Differential cross-sections 
for inclusive 7r°-meson production are presented as a function of Bjorken-x and the four- 
momentum transfer Q 2 , and as a function of transverse momentum and pseudorapidity. A 
recent numerical BFKL calculation and predictions from QCD models based on DGLAP 
partem evolution are compared with the data. 
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1 Introduction 



Hadronic final state analyses in Deep-Inelastic Scattering (DIS) interactions at HERA allow 
novel stringent tests of the physics of Quantum Chromodynamics (QCD), the theory of the 
strong interactions, in a kinematical region of high parton densities which so far has not been 
accessible [jT], 0, |3|, |]. The high Center of Mass System (CMS) energy (~ 300 GeV) of the 
HERA collider allows a region in Bjorken-x of ~ 1CT 5 — 10~ 4 to be reached while keeping 
the momentum transfer, Q 2 , larger than a few GeV 2 , hence remaining in the regime of pertur- 
bative QCD (pQCD). In DIS a parton in the proton can induce a QCD cascade consisting of 
several subsequent parton emissions before the final parton interacts with the virtual photon. 
The multiplicity and the x distribution of these emitted partons differ significantly in different 
approximations of QCD dynamics at small x. 

At low x, pQCD evolution is complicated by the occurrence of two large logarithms in the 
evolution equations, namely In 1/x and In Q 2 . In contrast, in the better tested region of pQCD 
at larger x a summation of the leading In Q 2 terms is sufficient. A complete perturbative treat- 
ment in the low-s region is not yet available, and different approximations are made resulting 
in different parton dynamics. At high Q 2 and high x pQCD requires the resummation of con- 
tributions of a s \n(Q 2 /Ql) terms, yielding the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi) [|5p evolution equations. However at small x the contribution of large leading In 1/x 
terms may become important. Resummation of these terms leads to the BFKL (Balitsky-Fadin- 
Kuraev-Lipatov) [|p evolution equation. Hence a pertinent and exciting question is whether 
these In 1/x contributions to the parton evolution can be observed experimentally. 

Differences between different dynamical assumptions for the parton cascade are expected 
to be most prominent in the phase space region towards the proton remnant direction, i.e. away 
from the scattered quark. Here we investigate the region of central rapidity in the CMS system 
of the virtual photon and the proton. In the HERA laboratory frame this corresponds to a region 
of small polar angles and has been generically termed "forward region"Q. In previous HI anal- 
yses results have been presented on forward jet and forward inclusive charged and neutral pion 
production [jT], ^J], based on data collected in 1994. A measurement on forward jet production 
has been presented by the ZEUS collaboration [Q. In this paper we study forward single %° 
production for a considerably larger data sample than that of [jSp, collected in 1996 and which 
allows the selection of particles with larger transverse momentum, px- The production of high 
Pt particles is strongly correlated with the emission of hard partons in QCD and is therefore 



sensitive to the dynamics of the strong interaction [ |T0| , [IT], [12]]. An advantage of studying single 
particles, as opposed to jets, is that no jet algorithm is needed and the potential to reach smaller 
angles than is possible with jets with broad spatial extent. Furthermore, theoretical calculations 
at the parton level can be convoluted with known fragmentation functions [JT^], allowing a direct 
comparison of the measurements and theory. The selection of high pr particles is also inspired 
by the proposal of Mueller [[]]] to select events where the photon virtuality Q 2 and transverse 
momentum squared of the parton emitted in the parton cascade, k\, are of similar magnitude, 
thereby suppressing the k T ordered DGLAP evolution with respect to the non-k T ordered BFKL 
evolution. 



1 HI uses a right-handed coordinate system with the z-axis defined by the incident proton beam and the y-axis 
pointing upward. 
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In this analysis 7r°'s are selected in DIS events at low x in the region of momentum transfer 
2 < Q 2 < 70 GeV 2 . The 7r°'s are required to have a polar angle in the lab frame between 
5° and 25°, and transverse momentum larger than 2.5 GeV in the hadronic CMS (contrary 
to the analysis in [§], where a minimum transverse momentum of 1 GeV was required in the 
laboratory frame). Large transverse momenta in the hadronic CMS, as opposed to the laboratory 
system, are more directly related to hard subprocesses, since in the quark parton model picture 
the current quark has zero in the hadronic CMS. The increased transverse momentum cut 
enhances the sensitivity to hard parton emission in the QCD cascade and provides a hard scale 
for perturbative calculations. It also reduces significantly the influence of soft hadronization. 

A calculation based on pQCD which uses the BFKL formalism for the perturbative part, 
and fragmentation functions for the hadronization, is available [|J and will be compared with 
the data. In addition, models using 0(a s ) QCD matrix elements and parton cascades according 
to DGLAP evolution, and colour string hadronization, will be compared with the data. 



2 Experimental Apparatus 

A detailed description of the HI detector can be found elsewhere [[14]]. The following section 
briefly describes the components of the detector relevant for this analysis. 

The hadronic energy flow and the scattered electron are measured with a liquid argon (LAr) 



calorimeter and a backward SPACAL calorimeter, respectively. The LAr calorimeter [ |15| ] ex- 
tends over the polar angle range 4° < 9 < 154° with full azimuthal coverage. It consists 
of an electromagnetic section with lead absorbers and a hadronic section with steel absorbers. 
With about 44 000 cells in total, both sections are highly segmented in the transverse and the 
longitudinal direction, in particular in the forward region of the detector. The total depth of 
both sections varies between 4.5 and 8 interaction lengths in the region 4° < 9 < 128°. 
Test beam measurements of the LAr calorimeter modules showed an energy resolution of 



a E /E « 0.50/ y/E [GeV] © 0.02 for charged pions and of a E /E « 0.12/ y/E [GeV] © 0.01 
for electrons [JT3|] . The hadronic energy measurement is performed by applying a weighting 
technique in order to account for the non-compensating nature of the calorimeter. The absolute 
scale of the hadronic energy is presently known to 4%. The scale uncertainty for electromag- 
netic energies is 3% [|16p for the forward region relevant for this analysis. 

The SPACAL JT7| ] is a lead/scintillating fibre calorimeter which covers the region 153° < 
9 < 177.8° with an electromagnetic section and a hadronic section. The energy resolution for 
electrons is 7.5 %/yE © 2.5 %, the energy resolution for hadrons is ~ 30%. The energy scale 
uncertainties are 1% and 7% for the electrons and hadrons respectively. The timing resolution 
of better than 1 ns in both sections of the SPACAL is exploited to form a trigger decision and 
reject background. 

The calorimeters are surrounded by a superconducting solenoid providing a uniform mag- 
netic field of 1. 15 T parallel to the beam axis in the tracking region. Charged particle tracks are 
measured in the central tracker (CT) covering the polar angular range 25° < 9 < 155° and the 
forward tracking (FT) system, covering the polar angular range 5° < 9 < 25°. The CT consists 
of inner and outer cylindrical jet chambers, 2-drift chambers and proportional chambers. The 
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jet chambers, mounted concentrically around the beam line, provide up to 65 space points in 
the radial plane for tracks with sufficiently large transverse momentum. 

A backward drift chamber (BDC) in front of the SPACAL with an angular acceptance of 
151° < 9 < 177.5° serves to identify electron candidates and to precisely measure their di- 
rection. Using information from the BDC, the SPACAL and the reconstructed event vertex the 
polar angle of the scattered electron is known to about 0.7 mrad. 

The luminosity is measured using the reaction ep — > ep'j with two TICl/TlBr crystal calorime- 
ters installed in the HERA tunnel. The electron tagger is located at z = —33 m and the photon 
tagger at z = — 103 m from the interaction point in the direction of the outgoing electron beam. 



3 Theoretical Predictions 

Predictions for final state observables are available from Monte Carlo models using 0(a s ) 
matrix elements and parton cascades according to the DGLAP evolution, and from numerical 
calculations based upon the BFKL formalism. In the following we describe the models and 
calculations used. 



3.1 Phenomenological QCD Models 

Implementations of 0(a s ) matrix elements complemented by parton showers based on the 
DGLAP splitting functions are available in the programs LEPT06.5 [[[§ and HERWIG5.9 
Jl9tl . The factorization and renormalization scales are set to Q 2 . The predictions of these mod- 
els should be valid in the region: a s (Q 2 ) \n(Q 2 /Q 2 ) ^ i an d a s (Q 2 ) ln(l/x) < 1. In LEPTO 



the Lund string model as implemented in JETSET7.4 [fZOJ is used to describe hadronization 
processes. LEPTO includes soft colour interactions in the final state which can lead to events 
with a large rapidity gap. HERWIG differs from LEPTO in that it also considers interference 
effects due to colour coherence and uses the cluster fragmentation model for hadronization. The 
versions of LEPT06.5 and HERWIG5.9 used consider only DIS processes in which the virtual 
photon is point-like. 

Recently a model has been proposed (RAPGAP2.06 [|2T|]) which is also based on the DGLAP 
formalism but includes contributions from processes in which the virtual photon entering the 
scattering process can be resolved. The relative contribution from resolved photon processes 
depends on the scale at which the virtual photon is probed. As in p2fl the factorization and 



renormalization scale in this paper is taken to be Q 2 + p\ {p\ of the partons from the hard 
subprocess). 



The model calculations in this paper were made with the CTEQ4M [E3p parton densities 



for the proton and the SAS-1D [24] parton densities for the virtual photon. QED corrections 



are determined with the Monte Carlo program DJANG06.2 Q25[]. The contribution of photons 
emitted in the forward direction from QED processes originating from the quarks in the proton, 
were found to be negligible. 
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In a previous paper [|]] we compared the results with ARIADNE [ |26| 1 and LDCMC pTj]. 
ARIADNE provides an implementation of the Colour Dipole Model (CDM) of a chain of in- 
dependently radiating dipoles formed by emitted gluons [[28]] . Unlike LEPTO, the cascade of 
the CDM is not ordered in transverse momentum. For the present analysis it was confirmed 
that the predictions depend strongly on the parameters controlling the "size" of the diquark and 
photon, and are therefore not explicitly compared with data in this paperQ. The linked dipole 
chain (LDC) model [E9T] is a reformulation of the CCFM [BOH equation, which forms a bridge 
between the BFKL and DGLAP approaches. Calculations of the hadronic final state based on 
this approach are available with the LDCMC 1 .0 Monte Carlo which matches exact first order 
matrix elements with the LDC-prescribed initial and final state parton emissions. The model 
however failed to describe the data in []8p and is therefore not considered further in this paper. 



3.2 BFKL Calculation 

Recently n° cross-sections have been calculated [Q] based on a modified BFKL evolution equa- 
tion in order 0(a s ) convoluted with n° fragmentation functions. The modified evolution equa- 
tions include the so called "consistency constraint" pT| , which limits the gluon emission at 
each vertex in the cascade to the kinematically allowed region. It is argued that this constraint 
embodies a major part of the non-leading m(l/x) contributions to the BFKL equation, which 
have been found to be very important The predictions of this modified BFKL equation are 
therefore expected to be more reliable than those without this constraint. The parton densities 
and fragmentation functions used in the calculation are taken from [Q] and [JT3|] respectively. 



In this paper we compare "set (iii)" of [Q] to the data. In this set the scale for the strong coupling 
constant as is taken to be the transverse momentum squared of the emitted partons, k\, and the 
infrared cut-off in the modified BFKL equation is taken to be 0.5 GeV 2 . Calculations with these 
parameters give a fair description of the forward jet cross-sections from [^p when taking into 
account hadronization corrections. The predictions are labelled "mod LO BFKL" in the figures. 



4 Measurement 



4.1 Data Selection 

The analysis is based on data representing an integrated luminosity of C = 5.8 pb^ 1 taken 
by HI during 1996. Deep-inelastic scattering events are selected and the event kinematics are 
calculated from the polar angle and the energy of the scattered positron. The four momen- 
tum transfer squared, Q 2 , and the inelasticity, y, are related to these quantities (neglecting the 
positron mass) by Q 2 = 4 E e E\ cos 2 y and y = 1 — ^ sin 2 y 1 , where E\ and E e are the 
energies of the incoming and the scattered positron respectively, and 9 e is the polar angle of the 
scattered positron. Bjorken-x is then given by x = Q 2 / (y ■ s), where s is the square of the ep 
center of mass energy. 

2 A good description of the data presented in this paper can be achieved (using ARIADNE4.10) e.g. choosing 
PARA(10)=1.7 and PARA(14)=1.0 [Hi. 
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Experimentally the scattered positron is defined to be the highest energy cluster, i.e. local- 
ized energy deposit, in the SPACAL with a cluster radius of less than 3.5 cm and an associated 
track in the BDC. Experimental requirements based on the energy and the polar angle of the 
scattered positron are used during the preselection but these are superseded by stronger kine- 
matic cuts which restrict the data to the range 0.1 < y < 0.6 and 2 < Q 2 < 70 GeV 2 . The 
restricted y-range ensures that the particles from the current quark are detected in the central de- 
tector, and not in the forward region, and that the DIS kinematics can be well determined from 
the measurement of the scattered positron. Photoproduction background is further reduced to 
a negligible level by requiring 35 < ^ ■ (Ej — p z j) < 70 GeV [ ]3~5] ] with Ej andp 2J the en- 



ergy and longitudinal momentum of a particle respectively, and where the sum extends over all 
detected particles in the event, except for those in the small angle electron and photon tagger. 
The reconstructed primary event vertex must have a z coordinate not more than 35 cm away 
from the nominal interaction point. The trigger is based on energy depositions in the SPACAL 
and demands multiple track activity in the central tracker. For the events used in this analysis 
the efficiency of this trigger is around 80%, determined using data from an independent second 
trigger. 

After the selection about 600K events are available for further analysis. 



4.2 Forward 7r° -Meson Selection 

A measurement of particle production at mid-rapidity in the hadronic CMS system requires 
small forward angles in the lab system. It is difficult to identify individual charged particles 
in the forward direction in an environment with a high density of charged particles. However 
the finely segmented HI LAr calorimeter allows the measurement of 7r°'s down to very small 
angles. They are measured using the dominant decay channel n° — > 27. The tt° candidates 
are selected in the region 5° < 9^ < 25°, where 9 n is the polar angle of the produced n° . 
Candidates are required to have an energy such that x w = E w /E prot on > 0.01, with E proton the 
proton beam energy (820 GeV), and a transverse momentum in the hadronic CMS, w , greater 
than 2.5 GeV. At the high n° energies considered here, the two photons from the decay cannot 
be separated, but appear as one object (cluster) in the calorimetric response. Therefore, the 
standard method to identify 7r°-mesons by reconstructing the invariant mass from the separate 
measurement of the two decay photons is not applicable. 

In this paper, a detailed analysis of the longitudinal and transverse shape of the energy depo- 
sitions is performed to separate electromagnetic from hadronic showers. This approach is based 
on the compact nature of electromagnetic showers as opposed to showers of hadronic origin, 
which are broader. The analysis of shower profiles is made possible by the fine granularity of 
the calorimeter in the forward direction. It has a typical lateral cell size of 3.5 x 3.5 cm 2 . This 
can be compared to the mean Moliere radius R m which is 3.6 cm and the mean radiation length 
X which is 1.6 cm. The calorimeter has a four-fold longitudinal segmentation for the electro- 
magnetic section which has a thickness of 20 to 25 radiation lengths Xq. The main experimental 
challenge in this analysis is the high activity in this region of phase space, with hadronic show- 
ers "masking" the clear electromagnetic signature from the n° — > 27 decay. The overlap of a tt° 
induced cluster with another hadron is mainly responsible for losses of ir° detection efficiency, 
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since the distortion of the shower shape estimators it causes will, in many cases, lead to the 
rejection of the cluster candidate. 

The reconstruction of LAr data is optimized to contain all the energy of an electromag- 
netic shower in one cluster JT4]]. A 7r°-meson candidate is required to be a cluster with more 
than 90 % of its energy deposited in the electromagnetic part of the LAr calorimeter. A "hot" 
core consisting of the most energetic group of contiguous electromagnetic calorimeter cells of a 
cluster, which must include the hottest cell, is defined for each candidate [^j. More than 50 % 
of the cluster energy is required to be deposited in this core. The lateral spread of the shower 
is quantified in terms of lateral shower moments calculated relative to the shower's principal 



axis Q36| ] and required to be less than 4 cm |]T2|] . The longitudinal shower shape is used as a 



selection criterion via the fraction of the shower's energy deposited in each layer of cells in the 
electromagnetic part of the calorimeter. The precise specifications of these layers can be found 
in [[^jp. The part of the cluster's energy measured in the second layer minus that measured in 
the fourth layer is required to be more than 40% of the total cluster energy. This selects showers 
which start to develop close to the calorimeter surface and are well contained in the electromag- 
netic part of the calorimeter, as expected for showers of electromagnetic origin. As mentioned 
above, with this selection one cannot distinguish photons from a 7r°-meson decay and photons 
from other sources. The high energy required in the selection, however, ensures that contribu- 
tions from sources other than high energy 7r°-mesons (such as prompt photon production) are 
at a negligible level [T37p. The influence of ?y-meson production is corrected for in the analysis. 



Uncertainties in the relative rj and tt° production rates in the Monte Carlo models used have 
been studied and were found to have an negligible effect on the results. 

With this selection about 1700 (600) ir° candidates are found in the kinematic range 5° < 
6 n < 25°, x v > 0.01 and p Tv > 2.5 (3.5) GeV, with a detection efficiency better than 45%. 
Monte Carlo studies, using a detailed simulation of the HI detector for a sample of DIS events, 
yield a purity of about 70% for the selected 7r°-meson sample. The impurities are due to 
misidentified hadrons and from secondary interactions of charged hadrons with passive ma- 
terial in the detector (between one and two radiation lengths in the forward region). These 
studies show that less than 10% of the selected n° candidates stem from secondary scattering 
of charged hadrons with passive material in the forward region, where the amount of material 
between the interaction point and the calorimeter surface is largest. 

The determination of the n° acceptance and purity depends only on the particle density and 
energies in the forward calorimeter. To ensure that the different Monte Carlo models used are in 
reasonable agreement with the data in this respect the transverse energy flow, Et, around the n° 
candidate clusters was studied in detail. Both the Et flow and the Et spectra in the tail of the Et 
flow distributions are reasonably well described by the models used to determine the detector 
corrections [|T2"|] . Of the two models used, however, ARIADNE showed a higher particle density 



while LEPTO has a lower particle density than the data. Remaining differences of the detector 
corrections determined with the two Monte Carlo models are therefore used to estimate the 
systematic error. When the energy and transverse momentum requirements are lowered, the 
two photons from the n° decay become separable and a clear n° mass peak can be observed 
which is also well reproduced by the HI detector simulation. The same method of selecting 
7r°-mesons as outlined above was used in a previous HI analysis [§], where a measurement of 
charged particles in the same region was also performed. The measured tt° cross-sections were 
found to agree well with the average of the n + and n~ cross sections. Furthermore, the results 
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of the analysis of the 1994 data are found to be in good agreement with the present analysis in 
their overlapping phase space regions. 



5 Results 



The experimental results of the analysis are presented as differential ep cross-sections of for- 
ward 7r°-meson production as a function of Q 2 , and as a function of x, r] w and p^ n in three 
regions of Q 2 for > 2.5 GeV. The pseudorapidity 77^ is given by — In [tan (9/2)) with 9 
being the polar angle of the n° in the laboratory frame. In addition the ir° cross-sections as a 
function of x and Q 2 are measured for data with the threshold of the ir° transverse momentum 
increased to p?^ > 3.5 GeV. An increased p* T w threshold is expected to enhance the sensitivity 
to hard parton emission in the parton cascade. The phase space is given by 0.1 < y < 0.6, 2 < 
Q 2 < 70 GeV 2 , 5° < 6 W < 25° and x w = E^jE^on > 0.01, in addition to the p?^ thresholds 
given above. 6 n , E w and E proton are measured in the HI laboratory frame; ptp n is calculated 
in the hadronic CMS. The measurement extends down to x = 4-10~ 5 , covering two orders of 
magnitude in Bjorken-x. 

All observables are corrected for detector effects and for the influence of QED radiation by 
a bin-by-bin unfolding procedure. The detector effects include the efficiency, purity and ac- 
ceptance of the 7r°-meson identification as well as contributions from secondary scattering in 
passive material. The correction functions are obtained with two different models (ARIADNE 
and LEPTO) and detailed detector simulation. The final correction is performed with the av- 
erage of the two models. The remaining background from photoproduction in the data sample 
has been studied using a sample of photoproduction Monte Carlo events (PHOJET [[381]) repre- 
senting an integrated luminosity of about 1 pb _1 . The contribution from such events is found to 
be negligible in all bins. 

The typical total systematic uncertainty is 15-25%, compared to a statistical uncertainty of 
about 10%. Contributions to the systematic error include: the uncertainty of the luminosity 
measurement (1.8%), the statistical uncertainty in the determination of the trigger efficiency 
(5%), the uncertainty of the electromagnetic energy scale of the LAr (3%) and the SPACAL 
(1%) calorimeters which each contribute 5-10%, the variation of 7r°-meson selection and accep- 
tance requirements within the resolution of the reconstructed quantities (5-10%), and the model 
dependence of the bin-by-bin correction procedure using differences between ARIADNE and 
LEPTO (5-10%). 

The cross-sections as a function of x, shown in Fig. [j](a), exhibit a strong rise towards small 
x. In this and the following figures the inner error bars give the statistical errors, while the outer 
error bars give the statistical and systematical error added quadratically. It is of interest to note 
that the rise in x in Fig. |T](a) is similar to the rise of the total inclusive cross-section as measured 
e.g. in [JJ^p. This is demonstrated in Fig. |T](b), which shows the rate of 7r°-meson production 



in DIS as a function of x obtained by dividing the cross-section shown in Fig. [T](a) by the 
inclusive ep cross-section in each bin of x and Q 2 . The inclusive cross-section is calculated by 



integrating the HI QCD fit to the 1996 structure function data as presented in [ p9| ] for every bin 
of inclusive 7r°-meson cross-sections. Note, however, that the n° rate increases with increasing 
Q 2 . The x-independence seen in Fig. |T](b) in a fixed Q 2 interval implies that the tc° rate for 
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particles with a px above the cut-off and within the selected kinematical region, is independent 
of W, the hadronic invariant mass of the photon-proton system. 

The shapes of da^/dr]^ and da n / dp T n (Fig. |2|) show no significant dependence on Q 2 . The 
measurements of the latter extend to values of transverse momenta as high as 8 GeV. Since r] n is 
measured in the laboratory frame, mid-rapidity in the hadronic CMS corresponds approximately 
to 7]^= 2 in Figure [Z] (a). Figure [3] shows the inclusive 7r°-meson cross-section as a function of 
Q 2 . The cross-section falls steeply with increasing Q 2 . Figure |] finally shows da n /dQ 2 and 
da n /dx for the higher threshold of p^ n > 3.5 GeV. No significant change in shape of the 
distributions occurs when the p^ n threshold is raised, but the cross-sections are reduced by 
about a factor of three. 

All differential cross-sections are compared to three predictions based on different QCD ap- 
proximations. The DGLAP prediction for pointlike virtual photon scattering (including parton 
showers) as given by LEPT06.5 falls clearly below the data. There is still a fair agreement with 
data for the highest x and Q 2 bins, shown in Fig. [I](a) and Fig. |3], but differences occur in the 
low-x region. These are as large as a factor of five in the lowest x region. LEPTO also fails to 
describe the ratio in Fig. [T](b), and shows a strong decrease with decreasing x. The mechanism 
of emitting partons according to the DGLAP splitting functions, combined with pointlike virtual 
photon scattering only, is clearly not supported by the data, in particular at low x. The LEPTO 
prediction is based on about seven times the integrated luminosity of the data. Comparisons to 
HERWIG5.9 (not shown) lead to similar conclusions [[l2|]. 



A considerable improvement of the description of the data is achieved by a model which 
considers additional processes where the virtual photon entering the scattering process is re- 
solved. This approach can be regarded as an effective resummation of higher order corrections. 
It also provides a smooth transition towards the limit of Q 2 = 0, i.e. photoproduction, where 
resolved processes dominate in the HERA regime. Such a prediction is provided by RAP- 
GAP2.06 [|2~T|]. In Fig. [T](a) RAPGAP2.06 predicts a cross-sections very close to the measured 
distributions with the exception of the lowest Q 2 bin where the prediction is too low. All pre- 
dicted cross-sections increase by up to 30% when the scale in the hard scattering is increased 
from Q 2 + p\ to Q 2 + [JL^], and therefore do not improve the overall description signif- 
icantly. Hence RAPGAP, with the parton distributions used here, does not describe the low-x 
behaviour of the data over the full range. The RAPGAP prediction is based on approximately 
four times the integrated luminosity of the data. 

The ARIADNE model (not shown), with parameters as given before, can describe the data 
presented in this paper [|T2|], but it remains to be shown whether this choice allows for a consis- 



tent description of other aspects of the DIS final state data. Moreover, a moderate variation of 
these parameters leads to large changes in the prediction. 

Next we compare the data with a prediction of the 7r°-meson cross-section based on a modi- 
fied LO BFKL parton calculation convoluted with n° fragmentation functions. The predictions 
obtained with these calculations turn out to be in good agreement with the neutral pion cross- 
sections measured in most of the available phase space, but are below the data at the lowest 
values of Q 2 . The calculation also describes well the ratio shown in Fig. |T](b), except possibly 
at the largest Q 2 , x bin. This ratio has been calculated by using the corresponding prediction for 



the inclusive cross section, based on the BFKL formalism [32]. 
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The BFKL predictions involve a cut-off parameter in the transverse momentum squared 
fc£ of the partons: k 2 =0.5 GeV 2 and a choice of renormalization scale. It is shown in [Q] that a 
variation of k$ by a factor of two leads to less than a 10% change of the cross-sections. The scale 
dependence is larger; a change from k\ tok^/A leads to an approximate increase of 60% of the 
cross-sections. This affects mostly the normalization, but not the shape of the distributions. 

The good agreement between this prediction and the data suggests that the modified BFKL 
evolution equation, using the consistency constraints, is a good approximation for low-x evo- 
lution in the considered phase space. This, in turn, can then be interpreted as a sign of the 
experimental manifestation of leading In 1/x terms which are anticipated in pQCD evolution. 

6 Conclusions 



Differential cross-sections of forward n° production have been measured for particles with 
Vt^ > 2.5 (3.5) GeV, 5° < 9 W < 25° and x v = E w /E proton > 0.01, for DIS events with 0.1 < y < 
0.6 and 2 < Q 2 < 70 GeV 2 . The data are sensitive to QCD parton dynamics at low x (high 
parton density) and mid-rapidity in the hadronic CMS system. They discriminate between 
different approximations to QCD evolution in the new regime opened up by HERA. The data 
show a strong rise of the cross section with decreasing x. This rise is similar to the rise of the 
inclusive cross section. 

Models using 0(a s ) QCD matrix elements and parton cascades according to the DGLAP 
splitting functions cannot describe the differential neutral pion cross-sections at low x. Inclusion 
of processes in which the virtual photon is resolved improves the agreement with the data, but 
does not provide a satisfactory description in the full x and Q 2 range. A calculation based on 
the BFKL formalism is in good agreement with the data, particularly for the shape description, 
but the absolute normalization remains strongly affected by the scale uncertainty. So far the 
data in the phase space selected in this analysis could not be confronted with a next-to-leading 
(NLO) order prediction - either for DGLAP or for BFKL. More definite conclusions therefore 
have to be delayed until such calculations become available. 
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Table 1: The numerical values of the inclusive differential 7r°-meson cross sections as shown in 
Figures 1-4. The value of the cross section is presented together with the respective statistical 
and total error. The phase space is given by 0.1 < y < 0.6, 5° < 9 n < 25° and x w > 0.01 in 
addition to the restrictions given in the Table. 



15 



-O 
c 



X 



1500 



cf 1000 
■o 

500 


400 
200 


100 
50 



H1 data 

P ; I[ >2.5 GeV 

■■ LEPT0 6.5 
- RAPGAP 2.06 

- mod. LO BFKL 



[2.0 ...4.5] 



I I I I I I I L 



Q 2 / GeV 2 
[4.5 ... 15.0] 



I I I I I I I L 



[ 15.0 ... 70.0] 



a) 



It 



i i i i 



DC 



0.003 
0.002 
0.001 


0.004 
0.003 
0.002 
0.001 


0.008 
0.006 
0.004 
0.002 





o R x (H1 do/dx) 
P ;, 7[ >2.5 GeV 

■■■■ LEPT0 6.5 

- - RAPGAP 2.06 

— mod. LO BFKL 

[2.0... 4.5] 

" I i i I i i 



[4.5 ... 15.0] 




I I I I I I l_L 



[ 15.0 ... 70.0] 

Q 2 / GeV 2 



b) 



I I I I I I l_L 



10 



10 



10 



-4 



10 



Figure 1: Inclusive 7r°-meson production cross-sections as a function of x for p* T n > 2.5 GeV 
in three regions of Q 2 (a). The phase space is given by 0.1 < y < 0.6, 5° < 9 n < 25° and x w 
= E^/Eproton > 0.01 and the Q 2 ranges given in the figure. The inner error bars are statistical, 
and the outer error bars give the statistical and systematical error added quadratically. 9 W and 
x w are measured in the HI laboratory frame, p* Tlx is calculated in the hadronic CMS. The QCD 
models RAPGAP (sum of direct and resolved contributions) and LEPTO are compared to the 
data. Also shown is the prediction of the modified LO BFKL calculation by Kwiecinski, Martin 
and Outhwaite. Figure (b) shows the rate of 7r°-meson production in DIS as a function of x 
obtained by dividing the cross-section shown in (a) by the inclusive ep cross-section in each bin 
of x and Q 2 . 



16 



-O 
CL 



-d 200 



Q 2 /GeV 2 



[2.0... 4.5] • H1 data 

P ;, 3t >2.5 GeV 



T3 



100 



-i 



J I I I I I I I I I I I I I I L 



200 



100 




[ 15.0 ... 70.0] 



200 



100 



a) 



LEPTO 6.5 
RAPGAP 2.06 

mod. LO BFKL 



J I I I I I I I I I I I I I I L 



1.5 



2.5 



3 



> 

CD 
O 

23 10 

CL 



-8- 10 



10 



10 




10 



10 



• H1 data 

p; n > 2.5 GeV 

± 



Q 2 / GeV 2 



[2.0 ...4.5] 

I I I I L 



[4.5... 15.0] ■ 

I I I I L 



=4- 



b) 



[ 15.0 ... 70.0] 

I I I I L 



4 5 6 7 8 



Figure 2: Inclusive 7r°-meson production cross-sections as a function of rj^ (a) and w (b) for 
Pt-k > 2.5 GeV in three regions of Q 2 . The phase space is given by 0.1 < y < 0.6, 5° < 
9 W < 25° and x w = E n / 'E proton > 0.01 and the Q 2 ranges given in the figure. The inner error 
bars are statistical, and the outer error bars give the statistical and systematical error added 
quadratically. 9 n and x w are measured in the HI laboratory frame, is calculated in the 
hadronic CMS. The QCD models RAPGAP (sum of direct and resolved contributions) and 
LEPTO are compared to the data. Also shown is the prediction of the modified LO BFKL 
calculation by Kwiecinski, Martin and Outhwaite. 
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Figure 3: Inclusive 7r°-meson production cross-sections as a function of Q 2 for n > 2.5 GeV. 
The phase space is given by 0.1 < y < 0.6, 5° < 9 W < 25° and x w = E^/Eproton > 0.01. The 
inner error bars are statistical, and the outer error bars give the statistical and systematical error 
added quadratically. 9 n and x n are measured in the HI laboratory frame, is calculated in 
the hadronic CMS. The QCD models RAPGAP (sum of direct and resolved contribution) and 
LEPTO are compared to the data. Also shown is the prediction of the modified LO BFKL 
calculation by Kwiecinski, Martin and Outhwaite. 
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Figure 4: Inclusive 7r°-meson production cross-sections as a function of x (a) and Q 2 (b) for 
Vt^ > 3.5 GeV. The phase space is given by 2.0 < Q 2 < 70.0 GeV 2 ,0.1 < y < 0.6, 5° < 
9 W < 25° and x n = E w / ' E proton > 0.01. The inner error bars are statistical, and the outer error 
bars give the statistical and systematical error added quadratically. 9 n and x n are measured in 
the HI laboratory frame, is calculated in the hadronic CMS. The QCD models RAPGAP 
(sum of direct and resolved contribution) and LEPTO are compared to the data. Also shown is 
the prediction of the LO BFKL calculation by Kwiecinski, Martin and Outhwaite. 
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